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ABSTRACT 

We introduce a graph-theoretic approach to synchro¬ 
nizing clocks in an ad hoc network of N timepieces. 
Clocks naturally drift away from being synchronized 
because of many physical factors. The manual way of 
clock synchronization suffers from an inherrent propaga¬ 
tion of the so called “clock drift” due to “word-of-mouth 
effect.” The current standard way of automated clock 
synchronization is either via radio band transmission 
of the global clock or via the software-based Network 
Time Protocol (NTP). Synchronization via radio band 
transmission suffers from the wave transmission delay, 
while the client-server-based NTP does not scale to 
increased number of clients as well as to unforeseen 
server overload conditions (e.g., flash crowd and time- 
of-day effects). Further, the trivial running time of 
NTP for synchronizing an iV-node network, where each 
node is a clock and the NTP server follows a single-port 
communication model, is 0{N). We introduce in this 
paper a 0(\ogN) time for synchronizing the clocks in 
exchange for an increase of 0(N) in space complexity, 
though through creative “tweaking,” we later reduced 
the space requirement to 0(1). Our graph-theoretic 
protocol assumes that the network is IK at, while the 
subset of clocks are in an embedded circulant graph 
^n<N with q jumps and clock information is commu¬ 
nicated through circular shifts within the All 

N nodes communicate via a single-port duplex channel 
model. Theoretically, this synchronization protocol 
allows for N{log N)~^ — 1 more synchronizations than 
the client-server-based one. Empirically through sta¬ 
tistically replicated multi-agent-based microsimulation 
runs, our protocol allows at most 80% of the clocks 
synchronized compared to the current protocol which 
only allows up to 30% after some steady-state time. 
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1. INTRODUCTION 

The “Juan Time, On Time” is a project of the Depart¬ 
ment of Science and Technology (DOST) launched in 
30 September 2011 which aim to campaign for the use 
of the Philippine Standard Time (PST). Since 1978, 
the PST is legally and officially maintained by DOST’s 
Philippine Atmospheric, Geophysical and Astronom¬ 
ical Services Administration (PAGASA) [3]. However, 
due to various reasons, the PST has not been utilized 
by Filipinos, whether in public or private transactions, 
resulting to having timepieces that are not synchronized 
with the PST. There are many problems that result 
by having non-PST-synchronized timepieces, some pos¬ 
sible (though relatively exaggerated) examples of these 
are: 


1. Historical and official events being recorded with 
conflicting times - e.g., in law enforcement, blot¬ 
ters with conflicting records of when crimes were 
committed may cause the criminal justice system 
to incarcerate an innocent person or free a guilty 
one. 

2. Financial transactions, specifically those done 
electronically, may cause one investment to lose a 
supposedly financial gain - e.g., an online bidder 
may submit a bid which might be a second late 
because heiQ timepiece is not synchronized with 
that of the bidding institution’s. 

3. In national defense, the order of a military com¬ 
mander may be executed several seconds earlier 
or later, instead of on time, which may later prove 
fatal to national security concerns - e.g., an air 
bomber pilot may release a second early a bomb 
payload to a rebel camp holding up hostages that 
have not yet been evacuated to a safe zone. 

4. In scientific research that rely on the accuracy 
and timeliness of the measuring devices - e.g, a 
clock-based data monitoring device may provide 
a sequence of wrong data array because the clock 
ran faster than expected, which if not corrected 
may prove crucial to the research conclusion. 


^Note: The use of the female gender in this paper is 
just a writing style and this could mean either without 
being prejudice to the other. 



1.1 Clock Drift 

There are many reasons why timepieces are not synchro¬ 
nized with one another, even though they started accu¬ 
rately ^nchronized. One of the reasons is the “clock 
drift” [1^ which happens because of the following phys¬ 
ical reasons: 

1. The clock changed its frequency (i.e., frequency 
shift). 

2. The clock changed its phase (i.e., phase shift). 

3. For a limited time (i.e., maybe a burst of sev¬ 
eral milliseconds), the clock experienced an 
unstable/interrupted power supply that resulted 
in either a frequency shift, a phase shift, or both. 

4. During an extended use or because of environ¬ 
mental factors, a clock was heated up that resulted 
in a frequency shift, a phase shift, or both. 

1.2 Word-of-mouth Propagation of Clock 
Drift 

In the past, and even until now, timepieces are generally 
updated using the following simple process: 

1. Query a supposedly trusted and authoritative time 
source, which usually is a person, a radio station 
announcing a time check, or a TV station showing 
time; and 

2. Manually reset the timepiece to the exact time 
returned by the time source, without considering 
the lag time between receiving the information 
from the source and the time it took to reset the 
clock. 

Because of this process, the recipient of the query 
answer would have reset her timepiece with an inherent 
“clock drift” due to “word-of-mouth” effect illustrated 
as follows: Given N persons pi,p 2 , • • • where pi is 
an official authoritative source of time. If p 2 updates 
her timepiece by querying pi , and then p 3 updates her 
timepiece by querying p 2 , and then so on in a linear 
fashion up to pw updating her timepiece by querying 
Pn- 1 , at the {N — l)th step, pN would have a clock 
drift with an optimistic factor of at most 2N. This 
factor is due to the “word-of-mouth” propagation of the 
time lag. 

1.3 Clock Synchronization via Radio Trans¬ 
mission 

In advanced countries where timeliness is of utmost 
importance, like the United States and Japan, time¬ 
pieces are equipped with (usually an amplitude modu¬ 
lation or AM) radio band receiver and are updated 
or synchronized at specified frequency by a signal from 
a dedicated (usually government-run AM) radio trans¬ 
mitter. The transmitter is connected to a time standard 
device, such as an atomic clock. Timepieces in these 
areas automatically adjust to differences in time zones, 
as well as to changes in daylight saving times (DST). 


However, timepieces are only adjusted up to a resolu¬ 
tion of a second, because the respective AM receivers 
are not equipped to detect for the propagation delay of 
the radio signal from the transmitters. On the average, 
the propagation delay is approximately 1 s for every 
300 Km distance the receiver is from the transmitter. 
Thus, this type of clock synchronization system is effec¬ 
tive only to timepieces that only require a resolution of 
up to a second, which currently are useful for general 
human use. 

1.4 Internet-based Clock Synchronization 
via the Network Time Protocol 

The Network Time Protocol (NTP) is a time syn¬ 
chronization protocol implemented in software for the 
purpose of synchronizing computer clocks over packet- 
switched, variable latency data networks, such as the 
Internet. The NTP uses a revised version of the 
Agreement Algorithm, also known as the Marzullo’s 
Algorithm 0 , to select time sources for estimating the 
accurate time from a poll of noisy sources. Time sources 
become noisy because of the effects of variable network 
latency, which the algorithm corrects by using a jitter 
buffer. The jitter buffer is computed earlier by profiling 
the round-trip times (RTT) of several zero-payload 
packets from a source node to a target node in the 
network. The time is synchronized via a hierarchical, 
semi-layered system of clock sources, starting from what 
is termed as Stratum 0, a device that is connected to 
an atomic clock. Stratum 1 devices are computers that 
are connected to Stratum 0 devices and normally act 
as servers for timing requests from Stratum 2 servers. 
In general. Stratum n devices connect to Stratum n — 1 
devices to synchronize time in a hierarchical client- 
server, master-slave fashion, where the masters are the 
devices in Stratum n — 1 and the slaves are the devices 
in Stratum n. In the Philippines, no Stratum 0 device 
has been officially established, even with the launching 
of DOST’s “Juan Time, On Time” campaign, which 
only uses the word-of-mouth propagation of the correct 
time with up to 1 minute resolution. Despite of this, 
most computer servers are potential Stratum 1 devices 
if they connect to known Stratum 0 devices abroad. 

1.5 Potential of Institutions as Statum 1 
Service Providers 

Nowadays, various local government and private institu¬ 
tions, particularly those in the highly urbanized areas, 
run several computer servers for providing ICT services 
to their constituents [T^. Some of these servers might 
be converted to run in dual-server modes with NTP. A 
dedicated cluster of NTP servers to act as a publicly- 
available Stratum 1 devices could be setup but may 
prove cost ineffective as more client computers connect 
and query the cluster for correct time at a higher res¬ 
olution and to synchronize clocks. With the expected 
improvement of telephone and communication services 
in the country [I|], particularly due to a healthier busi¬ 
ness competition that the ASEAN integration in 2015 
will bring UM, it is expected that the use of mobile 
computers among constituents will double every year. 
For a relatively small central business district with a 
pessimistic maximum estimate of 10,000 constituents. 


each owning at least one mobile computer that query the 
cluster for correct time, the cluster will be overwhelmed 
with answering queries for RTTs than for answering 
queries about the correct time. Thus, it is seen that the 
NTP is not an efficient protocol for synchronizing the 
devices beyond Stratum 1 for a very, very large client 
base. 


1.6 The Solution: Peer-to-peer Protocol 
for Synchronizing Clocks Beyond 
Stratum 1 

The problem with using NTP beyond Stratum 1 is that 
it uses a master-slave type of communication, where the 
master could be overwhelmed by slaves that number in 
tens of thousands, especially if the bandwidth does not 
scale with the increase of estimated users. With a con¬ 
stant bandwidth towards the master, it is necessary that 
the bandwidth used for answering RTTs and queries be 
distributed among the participating slaves via what is 
called a peer-to-peer (P2P) communication approach, 
similar to the strategy employed by the famous Bit Tor¬ 
rent protocol im. Thus, a new protocol is needed to 
query time and synchronize clocks for devices beyond 
the Stratum 1 device. 

We present in this paper an integrated knowledge in 
Process Theory and Graph Theory, particularly that of 
circular-shift process over circulant graphs C [IJ, to 
design a protocol for synchronizing N clocks in a com¬ 
plete network Kn and to show that the (logA/')-step 
protocol is correct and achievable. We show that our 
clock synchronization protocol is faster by a factor of 
logN", where N is the number of timepieces that are 
concurrently synchronizing. 

2. IMPROVED BERKELEY PROTOCOL 
WITH RECURSIVE DOUBLING TECH¬ 
NIQUE 

In the Berkeley Protocol (BP), given N clocks namely 
Co,Ci,... ,Cn-i with time readings To,Ti,...,T at-i, 
respectively, where To 7 ^ Ti ^ ^ Tn-i, the problem 

is to synchronize the times without relying on a global 
clock r. BP does this by averaging the N time read¬ 
ings with the assumption that no time reading is too 
extreme to effect a skew to the average. This can be 
performed in two ways, through an elected leader and 
through distributed computation. In the first method, 
an elected leader, usually Co, collects the respective 
N — 1 time readings, computes the average T, and then 
distributes T to N" — 1 others. In the second method, 
everybody broadcasts their own time readings to others, 
and they respectively compute the average without any¬ 
more additional communication. 

2.1 The Elected Leader Computes 

In the first method, the collection of the respective time 
readings takes N — 1 steps, as the leader Co needs to 
retrieve the time readings of Ci, C 2 ,..., and Cn-i one 
at a time. To compensate for the elapsed time due to 
collection of each time readings, every time a reading T 
is received, Co puts its own timestamp To,z on it. At 
the end of the N — 1 collection steps, Co would have 
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Figure 1: The timeline of Co showing when in 
Co’s own perspective of time it received the 
respective time readings, as well as when it com¬ 
puted T. 


collected the time readings Ti,T 2 ,..., and Ta^-i with 
respective timestamps To,i,To,2, • • •, and To,w-i. At 
the time of computation, which interestingly is at To, 
the zth time reading would have aged To — To,z, thus T 
must be corrected with this difference. Figure [T] shows 
the timeline of Co with respect to the receipt of the time 
readings at the respective timestamps. 

2.1.1 Computation of the T 
The average time T is computed as follows: 


T = 


T = 


T = 


In the above equations, it would have sufficed to stop 
with Equation [ 1 ] but we will soon see that the form in 
Equation [2] is practically useful in optimizing the space 
complexity of the methodology. The space complexity 
requirement of this method is discussed further below 
(Subsection 12.1.31) . 

It would have taken To,c time to compute for T, thus 
T must be corrected with this amount of computation 
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time also. After correction, T + To,c will be distributed 
by Co to the N — 1 other clocks. This will be done 
by Co one clock at a time for a total of — 1 steps, 
where each step, the elapsed time due to the previous 
communication will be added to the corrected T. Thus, 
Cl will receive T + To,c, C2 will receive T + To,c + T)o,i, 
where Oo,i is the elapsed time when Co sent the new 
time reading to Ci, C3 will receive T+To,c + T)o,i+Oo,2, 
where Oo,2 is the elapsed time when Co sent the new 
time reading to C2, and so on. In general, the zth clock 
will receive T + To,c + X]j =2 VI < z < iV. 

2.1.2 Time Complexity Requirement 

This method takes N —1 steps to collect the respective 
time readings, one step to compute for the average, and 
N — 1 steps to distribute the corrected average for a 
total of 2N — 2 steps. Thus the time complexity of this 
method is 0{N). 

2.1.3 Space Complexity Requirement 
Intuitively, Co needs N — 1 spaces to hold the N — 1 
collected time readings. This is what Equation [T] pro¬ 
vides at a glance. However, Co can just use 2 spaces to 
separately hold the running sum of the collected time 
readings and the running sum of the timestamps. This 
is what Equation!^ is showing. Co can reuse one of the 
two spaces to hold the corrected T. Thus, this method’s 
best space complexity is C(l). 

2.1.4 The Pitfall of Simplicity 

Regardless of the time and space complexities, the 
method suffers from simplicity because it did not con¬ 
sider the additional time it will take for the time 
readings to reach Co from their respective clocks. 
In Eigure \T\ above, Ti is basically the same as To,z, 
VO < z < A". This is not the case, however, because 
each clock either runs faster or slower than Co. When 
Co collects data from C^, it must have recorded the 
timestamp so,i at the start of its communication with 
Ci. Upon receipt of the time reading Ti from Ci, Co 
must have also recorded the timestamp To,i. If Ci is 
synchronized with Co, definitely so,z < Ti < To,i. If 
we assume that the time it takes for a request from 
Co to reach Ci is the same as the time it takes for 
the response from Ci to reach Co, then that time is 
Eo,z = 0.5(so,z + To,z)- The amount (so,i+To,i) is known 
in the literature as the roundtrip time (RTT) [2,0,El], 
and therefore RTTo,i = 2Eo,i. This amount is the one 
missing in the above discussion. Eigure [2] shows the 
visualization of these time values between the exchange 
of Co and Ci. 


Considering the asynchronous nature of the clocks, we 
can now obtain an estimate for Ti that is closer to its 
correct value and it is given as Ti + 0.5RTTo,z. With 
this corrected value. Equation must also be corrected 
into: 
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Figure 2: The respective timelines of Co and Ci 
showing the time values elapsed when initiating 
at timestamp so,z and completing at timestamp 
To,z the collection of Ti from C — i. 

2.1.5 Computation of the RTT 
Where will the RTT’s come from? Here, we propose a 
methodology that minimizes the error of the estimate 
for RTT. The reason for the error is that the time it 
will take for Co’s request to reach Ci is almost always 
not the same as the time it will take for CiS response 
to reach Co- Depending on the implementation of the 
communication protocol, Co’s initial request might as 
well go as little as one bit in length, say the value 0 
upon receipt by Ci to mean that the elected leader, 
Co, is requesting Ci to send its time reading Ti. The 
response, however, could involve a 32-bit integer, rep¬ 
resenting the number of seconds since some reference 
year. The propagation of a 1-bit data is faster than the 
propagation of a 32-bit data, especially to bandwidth 
constrained communication channels. Thus, we want to 
create a methodology that ensures that the RTT is rel¬ 
atively constant during the time of the collection of the 
time readings, and at the same time, we want to read 
Ti while this RTT is seemed to be non-changing. 

We propose the following algorithm: 

Algorithm 1: Computation of RTTo,z with Ti 


1. Set j = 0. 

2. Repeat the following: 

(a) Increment j by 1 

(b) Co sends a 1 to Ci at time so,i 

(c) Co receives a 32-bit long data from Ci at time 

Tb,i 

(d) Co computes for the RTTo,i,j = To,i — so,i 
until j = some statistically possible value 

3. Compute for the average RTT^^ = ^k=i ^^TTo, 
and its standard deviation afi. 

4. If (Jq i is within some set allowed threshold, then 
we move to step (5] else we go back to step \T\ 
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5. Co sends a 0 to Ci at time so,i 

6 . Co receives the 32-bit long Ti from Ci at time To,i 

7. We set j = 0 and repeat the steps in[^to collect 
j RTTo,z,j’s. 

8 . Compute for the average RTT^^ — j—^ ^^1=1 RTTo,z,a: 
and its standard deviation cr^^. 

9. If |RTT^^ — RTT^^I < some threshold and |cr^^ — 
cr^il < some threshold, 

• then Co accepts Ti with RTT^^, 

• else we repeat the whole process from step [U 

We want to set j in Algorithm 1 such that the time it 
takes to compute for the RTT will not dominate the 
time it takes to exchange the respective T’s. Unfor¬ 
tunately, j will depend on the state of the underlying 
network which can only be set through experience. We 
assume, however, that the network will not be a factor 
and that we can set j to a value that can provide a sta¬ 
tistically acceptable degree of freedom. We then further 
assume that the contribution of this algorithm to both 
the leader computes and the distributed computation 
approaches is C(l). 

2.1.6 Improvement of the Steps in Collecting 
Time Readings 

The collection of time readings in the original protocol, 
as shown in subsection l2.1.2l takes N—1 steps, or a time 
complexity of 0{N). We improved this time complexity 
to 0{logN) by utilizing a recursive doubling technique 
which we illustrate here with A" = 8 as follows. The 
procedure is completed in 3 steps, instead of seven steps. 

At step 1 , Co sends a 0 to C4. The 0 bit sent by Co will 
be propagated first to all clocks, while clocks which have 
already received the bit will participate in sending. At 
step 2 , Co sends a 0 to C2, while C4 propagates the 0 
to Ce. At step 3, Co sends a 0 to Ci, while at the same 
time C2 propagates the 0 to C3, C4 to C5, and Ce to 
C7. After step 3, all clocks would have received the 0 
from Co. 

The sending of the respective time readings will be done 
in the opposite manner, also in three steps as follows: 

At step 1 , Co receives Ti from Ci, and at the same time, 

C2 receives T3 from C3, C4 receives T5 from C5, and Ce 
receives Tr from C7. All pairs will follow the procedure 
outlined in Algorithm 1 . At step 2 , Co receives T 2 and 
the corrected T3 from C2, while C4 receives Tq and the 
corrected T 7 from Ce, again both utilizing Algorithm 1 . 

At step 3, using Algorithm 1 , Co receives T4, T5, Tq, 
and T7 from C4. 

In general, time readings are collected by Co via a 
recursive doubling method in C(log A) steps. However, 
the space complexity has increased to a corresponding 
(!)(logA) from C(l). Notice that the amount of data 
being transferred from Ci to Co doubles every step. 
Since the total number of doubling is also log A for 
A clocks, then the maximum amount of data to be 



Figure 3: An example circular 2-shift on a C|, 
which can be done via a series of two circular 
l-shift operations. 

passed is log A times of the original one. This max¬ 
imum happens in the last step, though. 

2.1.7 Distributing T in 0(\ogN) Time 
After Co has computed the T, it will distribute the 
average time to A — 1 clocks via the same recursive dou¬ 
bling technique. The corresponding time complexity is 
0(\ogN) while the space complexity is 0(1). 

2.2 Distributed Computation of T 

In the second method, each of the clocks Co, Ci, ..., and 
Cn-1 will collect time readings To,Ti,..., and Tn-i 
from the respective other clocks. Once this collection is 
completed, each of the clock will perform the averaging 
on their own, without any more further communication 
to the other clocks. Thus, our analysis focuses on a par¬ 
ticular distribution scheme for the time readings. Intu¬ 
itively, each clock can perform a collection of time read¬ 
ings from other clocks, one at a time. That is, each clock 
will be elected as a leader, collect the time readings, 
and then compute the average for itself without sharing 
it. If this is done in lexicographic way, and since we 
have already shown earlier in Subsection 12.1.21 that this 
particular method takes 0(N) time complexity, then 
this method, intuitively will cost 0(N log N) time. The 
question to be asked, then, is can we do better than 
this? 

2.2.1 The Circular Shift Operation 
Given a set of A nodes Ui, U 2 , • • •, Vn that form a 
regular circulant graph of order A with q jumps 
(or simply CU) [^, , the circular ^-shift opera¬ 
tion [ 13 , [H: is a special permutation of the nodes’ 

indexes such that node Vi sends a data packet to node 
V(i+g) mod N (Figure [3|). Researchers have long proved 
that the optimal number of steps for a circular g-shift 
on a is mm(q, N — q). To improve the performance 
of the distribution methodology discussed in Section 
above, we have to assume that the clocks are arranged 
in a Cjv- This is not impossible to do since any will 
perfectly embed into a Kw 0. 

Intuitively, the distribution of the time readings to all 
clocks only needs a circular (A — l)-shift operation, 
which only requires min (A — 1, A — (A — 1)) = 1 oper¬ 
ation. One can argue that this is true because a circular 















(N — 1)-shift operation is equivalently a circular (—1)- 
shift operation (i.e., a circular 1 -shift operation in the 
opposite direction). This is not the case, however, as we 
will soon see in our modification to the circular g-shift 
operation discussed in the next subsection. 

2.2.2 Circular (N - l)-Shift-Copy Operation 
We use the fact that a circular [N — 1)-shift operation 
can be done by a series of {N — 1) circular 1-shift opera¬ 
tions. We modified, however, each circular 1-shift oper¬ 
ation such that the receiving clock copies the time read¬ 
ings that has been shifted to it. We call our new opera¬ 
tion as a Circular g-Shift-Copy Operation. The circular 
{N — l)-shift-copy operation is simply a series of — 1 
alternating circular 1 -shift and copy operations. 

If a copy operation takes 1 step, then our circular g-shift- 
copy operation takes 2{N — 1) steps, or a complexity of 
0{N)^ a vast improvement to the intuitive time com¬ 
plexity discussed in Section [221 which is 0{N log N). 
Since each circular 1 -shift-copy operation only requires 
sending 1 data item per operation, then the circular 1 - 
shift-copy operation takes a space complexity of ( 2 ( 1 ). 
However, the receiving node must allocate a buffer that 
is equal to the amount of data that will be shifted, so 
the operation can take a maximum space complexity 
of 0{N). We can strategically reduce this space com¬ 
plexity to ( 2 ( 1 ) if for every intermediate circular 1 -shift- 
copy operation, the sum of the copied time readings will 
already be computed. 

Figure m shows a visualization of the progression of the 
first three 1-shift-copy operations on an iV = 8 clock 
network. 

2.2.3 Circular (N- l)-Shift-Copy Operation with 
Recursive Doubling 

The circular {N — l)-shift-copy vastly improves the time 
complexity of the operation from 0{N log N) down to 
0{N). The next question is, can we do better? It turns 
out that the answer to the question is a resounding yes 
as we shall soon see with our new proposed method we 
called recursively-doubled circular {N — 1)-shift-copy. 
This method takes the time complexity of (2(logiV) 
steps, which we will describe as follows: 

1. During the first step, instead of assuming that the 

clocks were arranged in a Cjv, we assumed that the 
clocks were arranged in a . This means that 

clock Ci will be connected to clocks Cij^yN/2\ and 
Ci- [Ar/2j, VO < z < A/". Such a circulant graph con¬ 
tains [A/'/2j disconnected C^’s. These subgraphs 
can alternately be seen as a linear graph L2 of 
order 2. The circular 1-shift-copy operation can 
be performed in these subgraphs concurrently. 

2. During the second step, we assumed that the 

clocks were arranged in a , where each 

clock Ci will be connected to clocks Ci^iN/4j 
and Ci-iN/ 4 }, VO < z < A^. Such a circulant 
graph contains [A^/4J disconnected Ci. As in the 
previous step, these subgraphs can concurrently 



Figure 4: An example progression of a circular 
7-shift-copy operation on a C|: (a) The data dis¬ 
tribution before the circular 7-shift-copy oper¬ 
ation; (b) The data distribution after the first 
l-shift-copy operation; (c) The data distribution 
after the second 1-shift-copy operation; and (d) 
The data distribution after the third 1-shift-copy 
operation. 


perform a circular 1-shift-copy operation each. In 
general, at step /c, we assumed that the clocks 
I N2~^ I 

were arranged in a Cjy ^ • The original network 
will be composed of [N2~^\ disconnected 
These subgraphs will concurrently perform a cir¬ 
cular 1-shift-copy operation each to distribute the 
data. 

3. At the last step (i.e., (logA)th step), the clock 
will be assumed to be arranged in a Cjv, where 
the circular 1-shift-copy operation distributes the 
final set of time readings. 


In this new method, the distribution of the time read¬ 
ings takes a time complexity of (2(log(A)). Figure [5] 
shows the evolution of the circulant graphs at each step 
of the methodology with N = 8. 


2.2.4 Complexities of the Recursively Doubled 

Circular (N - l)-Shift-Copy 
Intuitively, the time complexity of the recursively dou¬ 
bled circular [N — l)-shift-copy operation is 0(logN). 
Notice however that the space complexity doubles every 
step, with the (logA/')th step taking N spaces. Obvi¬ 
ously, the space complexity is (2(A). However, we can 
strategically reduce the space complexity if at every step 
we already compute the sum of the time readings. Thus, 
the space complexity of our proposed method can be as 
good as (2(1) space. 











Figure 5: An example progression of a recur¬ 
sively doubled circular 7-shift-copy operation on 
a C|: (a) The data distribution after step 1 where 
the clocks were arranged as a C|; and (b) The 
data distribution after step 2 where the clocks 
were arranged as a C|. 

2.3 Comparison between Leader Com¬ 
putes and Distributed Computation 

Both the Elected Leader Computes and the Distributed 
Computation have the same time and space complex¬ 
ities of 0(\ogN) and 0(1), respectively. This does 
not mean that we can now select any of the two in 
the implementation of the methodology. For all prac¬ 
tical purposes, we chose to implement the distributed 
computation method because the elected leader com¬ 
putes method will suffer from being “orphaned” when 
the elected leader decided to leave the network at the 
middle of the computation. Thus, we see the distributed 
computation method as a more robust method from the 
dynamism brought about by the constant movement of 
the clocks in and out of the ad hoc network. 


2.4 Implementation of the Distributed 
Computation through Computer Net¬ 
work Simulation 

We implemented the distributed computation by 
writing a program V that averages the internal clocks 
of computers connected in a local area network (LAN). 
We used a simple socket programming [2^ so that clock 
information can be distributed among the computers in 
the LAN, using the efficient recursively doubled circular 
{N — l)-shift-copy operation. 

Figure [6] shows the screen capture of an [N — 6)- 
clock synchronization problem implemented in six 
LAN-connected x686 processors, each running a multi¬ 
programming Gnu/Linux operating system. 

3. TIME SYNCHRONIZATION WITH 
MULTI-AGENT SYSTEM 

Our multi-agent-based [T^, [U time synchronization 
protocol is basically composed of three simple steps 
that can be implemented by any mobile simulated clock 
Ci depending on what state is in: (1) Update own 
time from the global clock T; (2) Update of own time 
from other clocks, and (3) Update other clocks. In our 
protocol, we assumed that each clock has the following: 


1. Time Record {Ti) - The current time of each clock; 
each clock has different time records since they are 
not synchronized; 

2 . r Synchronization Record (TSR^) - A time record 
when clock Ci last synchronized its time from T; 
and 

3. Clock Synchronization Record (CSR^) - A time 
record when Ci last synchronized its time from 
other clocks. 

The clock Ci can be in any of these two states: 

1. IN(r) - This means that Ci is under the influence 
of a global clock T; and 

2. OUT(r) - This means that Ci is not under the 
influence of T. 

3.1 Time Syuchrouizatiou uuder iN(r) 

When Ci is under the influence of a global clock T, 
Ci immediately synchronizes with T via a peer-to-peer 
protocol (P2P). The immediacy of the synchronization 
scheme assures those clocks which enter the circle of 
influence in an almost tangent to the edge of synchro¬ 
nization. Entering at a tangent means that these clocks 
will soon be out of the influence of T. Whenever Ci 
synchronizes with T, it updates with it its TSR as well. 
After the first synchronization, Ci may either be in 
one of the two available modes: (1) Passive Mode; or 
( 2 ) Aggressive Mode. These modes were developed to 
favor those clocks which are equipped with ranging- 
capable device. A clock with no ranging capability auto¬ 
matically chooses the aggressive mode, while a clock 
with ranging capability first chooses the passive mode 
and then switches to aggressive mode. When a range- 
capable clock Ci can range that its distance from T is 
decreasing, it uses the passive mode. However, when 
Ci can sense that its distance from T is increasing, then 
it switches to the aggressive mode. The passive mode 
allows for the conservation of power, especially for those 
clocks that are powered by batteries. 

3.1.1 Passive Mode at IN(r) State 

Upon entry of Ci into the influence of T, it first queries 
the r which always returns the current global time C. If 
\C — Ti\ is lesser than some threshold value Th^ then Ci 
does not do anything. However, the moment \C — Ti\ > 
Th, Ci immediately updates its Ti with G, as well as its 
TSR. While still under the influence of T, Ci continually 
queries the T for G, until |G — T^j > Th. 

3.1.2 Aggressive Mode at IN(r) State 
Regardless of the |G — T^j compared to T/i, Ci always 
immediately updates its Ti with G. The aggressive 
mode assures the clock that it always has the most 
recent TSR upon leaving the influence of T. 

3.2 Time Synchronization under ouT(r) 

When a clock Ci is out of the influence of T, then it 
could be under the influence of other clocks Gj, Vj ^ i 



Figure 6: Screen capture of executing the clock synchronization application V running on a 6-coniputer 
LAN through a remote secured shell (SSH) session. Shown in this screen capture are six terminals, 
each connected to the different computers where internal clock of each is being synchronized. 


within its immediate broadcast vicinity. Assuming that 
Ci enters a broadcast vicinity of iV — 1 other clocks, 
then we can use the protocols discussed in Section 
particularly the distributed computation scheme. How¬ 
ever, we will modify the protocol to compute for the 
max^Q^(rSRz) instead. We now propose a new method 
we called recursively doubling circular (A—l)-shift-max 
operation, where at each step of the operation, Ci com¬ 
pares its rSRi with what it received from its imme¬ 
diate neighbor, and retains the maximum between the 
two. This operation runs in 0(\ogN) time complexity 
and, since we only need to get the maximum, dehnitely 
with 0(1) space complexity. Figure [3 shows the visu¬ 
alization of the progression of time synchronization of 
N clocks using the recursively doubling (N — l)-shift- 
max method. 

For clocks with ranging capabilities, they will select to 
include those clocks that approach them into the net¬ 
work to prolong the life of their ad hoc community. Def¬ 
initely, those clocks that are already going away from 
them will soon be out of their group’s circle of influence. 
We do not want to include those clocks which may leave 
the network before the synchronization is completed. 


3.3 Simulation of the Protocol 

This protocol was simulated using a multi-agent-based 
simulation environment @5 05 [IS? [13 • We considered 
three scenarios as follows (Please refer to Figure [8|): 


1. Scenario A - In this scenario, we located the global 
clock (green circle) at the middle of the environ¬ 
ment, and placed three synchronization-disrupting 
areas (blue circles). Clocks are symbolized by the 
person icons, which randomly roam about the 
environment. 

2. Scenario B - This scenario is similar to Scenario A 
with the difference that the global clock is inside a 
fenced area and only those authorized persons are 
allowed to enter the area. This simulates the situ¬ 
ation wherein the global clock is only available to 
a few select people and that time synchronization 
will only happen if these select people will come 
in contact with those that were not selected. 

3. Scenario C - This scenario is similar to Scenario B 
but this time the fenced global clock is already 
located at the center, while the synchronization- 
disrupting areas are placed near the fence of the 
global clock. 

For comparison purposes for each scenario, we imple¬ 
mented a simple protocol that mimics how the current 
time-synchronization is currently being implemented. 
Synchronization happens when a newly F-synchronized 
clock Cx meets another clock Ci. In this protocol, 
Cx always shares its time with Ci via a simple P2P 
data exchange. Figure [9] shows the percentage of 
F-synchronized clocks within the hrst 30-s of the simu¬ 
lation time. This figure shows that the synchronization 
protocol that we developed can provide about 70% to 
































Figure 7: An example progression of a recur¬ 
sively doubled circular 7-shift-max operation on 
a C|: (a) The data distribution before step 1; 
The respective data distributions after steps 1 
(b), 2 (c), and 3 (d). 
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Figure 8: Snap shots of the multi-agent imple¬ 
mentation of the second protocol using a sim¬ 
ulation environment: (a) Scenario A; (b) Sce¬ 
nario B; and (3) Scenario C. 



Figure 9: Plot of the percentage of F- 

synchronized clocks during the first 30-s of the 
simulation for Scenario A (top line plot), Sce¬ 
nario B (middle line plot), and Scenario C 
(bottom line plot). The horizontal axis is 
in seconds while the vertical axis is in per¬ 
centage of F-synchronized clocks. Blue lines with 
square points are for the proposed protocol while 
orange lines with diamond points are for the 
simple protocol. 


80% synchronous clocks while the simple protocol can 
only provide up to 30% synchronous clocks for any 


scenario. 






















4. CONCLUSION 

In this paper, we argued that the DOST’s “Juan Time, 
On Time” program of using the PST with a simple syn¬ 
chronization protocol does not provide high percentage 
of T-synchronized clocks because of the inherrent clock 
drift brought about by the simple protocol. In fact, 
the clock drift is even enhanced by the simple protocol. 
We then provide an alternative automated protocol that 
synchronizes N clocks in 0(\ogN) time using only 0(1) 
memory. To prove that the proposed 0{\ogN) pro¬ 
tocol can provide a higher percentage of T-synchronized 
clocks, we simulated three scenarios where the proposed 
protocol is used. We compared the percentage of T- 
synchronized clocks to the same scenarios but this time 
when the simple protocol is used. For all scenarios, the 
proposed protocol provides 70% to 80% T-synchronized 
clocks while the simple protocol can only provide 20% 
to 30% T-synchronized clocks. Our protocol improved 
the number of T-synchronized clocks by at most 400% 
during the same time span. 
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